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Abstract 

We investigate whether Gaia can specify the binary fractions of massive stellar populations in the 
Galactic disk through astrometric microlensing. Furthermore, we study if some information about 
their mass distributions can be inferred via this method. In this regard, we simulate the binary astro- 
metric microlensing events due to massive stellar populations according to the Gaia observing strategy 
by considering (i) stellar-mass black holes, (ii) neutron stars, (iii) white dwarfs and (iv) main-sequence 
stars as microlenses. The Gaia efficiency for detecting the binary signatures in binary astrometric 
microlensing events is ^ 10 — 20 per cent. By calculating the optical depth due to the mentioned 
stellar populations, the number of the binary astrometric microlensing events being observed with 
Gaia with detectable binary signatures, for the binary fraction about 0.1, is estimated as 6, 11, 77 and 
1316 respectively. Gonsequently, Gaia can potentially specify the binary fractions of these massive 
stellar populations. However, the binary fraction of black holes measured with this method has the 
large uncertainty owing to a low number of the estimated events. Knowing the binary fractions in 
massive stellar populations helps for studying the gravitational waves. Moreover, we investigate the 
number of massive microlenses which Gaia specifies their masses through astrometric microlensing of 
single lenses toward the Galactic bulge. The resulted efficiencies of measuring the mass of mentioned 
populations are 9.8, 2.9, 1.2 and 0.8 per cent respectively. The number of their astrometric microlens¬ 
ing events being observed in the Gaia era in which the lens mass can be inferred with the relative error 
less than 0.5 toward the Galactic bulge is estimated as 45, 34, 76 and 786 respectively. Hence, Gaia 
potentially gives us some information about the mass distribution of these massive stellar populations. 

Keywords: Gravitational leasing: micro - Astrometry- Stars: binaries: general- methods: numerical 


1. INTRODUCTION 

Gaia is a space observatory of the European Space 
Agency (ESA), performing high precision astrometry, 
multi-colour, multi-epoch photometry and spectroscopy 
on astronomical objects brighter than 20 G magnitudes. 
This satellite, surveying the throughout sky, will observe 
each obj ect 72 times aver agely during five years of its 
lifetime (|Ever et al. 2013ll . The astrometric precision 
of Gaia depends strongly on the source magnitude 
and will be 30/ras to 600/ias for a single astromet- 
ric measurement o f stars with magnitude 13 to 19 
(|Varadi et al. 200^ . Different astronomical phenom¬ 
ena could b e observed during the G aia mission, e.g. 
supernovae (jBelokurov fc Evans 20031) , novae, dwarf 
novae, eruptive s tars, variable stars (|Varadi et al. 20091 
[Ever et al. 2013f), astrometric microlens ing events 
(jBelokurov fc Evans 2002t iProft et al. 201lh . etc. Here, 
we study some features of detecting the last one with 
Gaia. 

If the light of a background star passes through the 
gravitation field of a collinear foreground object, accord¬ 
ing to the Einstein general relati vity, the light beam 
bends toward the center of gravity (|Einstein 1936D . Con¬ 
sequently, several distorted, (un)magnified images are 
produced whose angular separations are too small to be 
resolved even by modern telescopes. Instead, the com- 
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bined and magnified light of images is received by ob¬ 
server, so-called gravitational microlensing. One of its 
feature is the displacement between the light centroid 
of images and the source position while the source star 
is passing the gravitational fi eld of lens, i.e. astromet¬ 
ric trajectory ( Walker 1995t iM ivamoto &: Yoshii 19951 : 
iHog et al. 1995H.leong et al. 19991) . During a microlens- 
ing event by measuring the astrometric leasing and par- 
allax effect, the mass of the deflector can be inferred 
(jPaczynski 19981 iMiralda-Escude 199611 . In microlensing 
events, the astrometri c cross-section is much larger than 
the photometric one (jPaczvhski 1996D . Hence, the op¬ 
tical depth of astrometric microlensing events is larger 
than the photometric one (jPominik fc Sahu 200(llj . 

Detecting the astrometric microlensing with the Gaia 
satellite was first discussed by Dominik & Sahu (2000). 
Then, Belokurov and Evans (2002) studied this sub¬ 
ject by more details.The astrometric microlensing due 
to massive objects in solar orbit beyond the Kuiper Belt 
was investigated by Gaudi and Bloom (2005). Also, Proft 
et al. (2011) investigated several catalogs to find several 
nearby stars with large proper motions which are poten¬ 
tial candidates for the astrometric microlensing detection 
in the Gaia era. 

Here, we investigate quantitatively the Gaia efficiency 
for specifying the hinarity fractions of microlenses in bi¬ 
nary astrometric microlensing due to massive stellar pop¬ 
ulations, i.e. stellar-mass black holes, neutron stars and 
white dwarfs in addition to main-sequence stars. In or¬ 
der to obtain the Gaia efficiency, we preform a Monte 
Garlo simulation. Then, we evaluate the optical depth 
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and the number of binary events being observed with 
Gaia with detectable binary signatures, considering dif¬ 
ferent binary fractions for massive stellar populations. 
Indicating the binary fractions in massive stellar popu¬ 
lations is useful for studying the gravitational waves and 
estimating their amounts, see e.g. Riles (2013), Cutler 
and Thorne (2002). We find that Gaia can potentially 
point out the binary fractions of massive stellar popu¬ 
lations. Furthermore, we estimate the Gaia efficiency 
for measuring the mass of massive stellar populations in 
the Galactic disk through astrometric microlensing of a 
single lens to indicate if Gaia can give any information 
about their mass distributions. Indeed, astrometric mi¬ 
crolensing events due to more massive microlenses have 
longer Einstein crossing times and larger angular Ein¬ 
stein radii. Therefore, the Gaia efficiency for measuring 
the lens mass in these events is high. 

This paper is organized as follows: In section m 
the basics of astrometric microlensing are explained. 
In the next section, we simulate binary astrometric 
microlensing events being observed with Gaia to in¬ 
vestigate how much information Gaia gives us about 
the binary stellar populations in the Galactic disk. In 
section o we study if Gaia can determine the mass 
distributions of massive stellar populations through 
detecting astrometric microlensing of a single lens. We 
conclude in the last section. 


2. BASICS OF ASTROMETRIC MICROLENSING 


In microlensing events, the light centroid vector of 
source star images does not coincide with the source 
position. This astrometric shift in source star posi¬ 
tion changes with time as a microlensing event pro¬ 
gresses. Eor a point-mass lens, the centr oid shift vec¬ 
tor of source star images is given by jWalker 19951 : 
iMivamoto fc Yoshii 1995t iHog et al. 19951 1: 


S0C 


fii 9 i + ^1262 
Ml + M2 


vlOe 


Oe 


( 1 ) 


where 9i and Mi are the position and magnification factor 
of fth image, u = pi -|- uqt) is the vector of the projected 
angular position of the source star with respect to the 
lens normalized by the angular Einstein radius of the lens 
in which p = {t — to)/I e , to is the time of the closest 
approach, is the Einstein crossing time, i and y are 
the unit vectors in the directions parallel with and normal 
to the direction of the lens-source transverse motion. The 
angular Einstein radius of lens is given by: 


Oe = v'» 

where Mi is the lens mass, k = AG/{c^A.U.) and Tr^ei = 
A.t/.(-^ — -j^) where Di and Ds are the lens and source 
distances from the observer. Usually, a microlensing par¬ 
allax is defined as tte = Ej-ei/dE which can be measured 
from the photometric event. According to equation 0, 
by measuring the relative parallax and the angular 
Einstein radius from astrometric observations the lens 
mass can be inferred. Nearby lenses have considerable 
and measurable relative parallaxes and the angular Ein¬ 
stein radii which are suitable candidates for astrometric 


measurements with Gaia. In that case, if the lens mass is 
high, the error bars of the lenses mass strongly decrease 
due to their large angular Einstein radii. 

In astrometric microlensing, the threshold amount of 
impact parameter i.e. Ua which gives an astrometric cen¬ 
troid shift more than a threshold amount i.e. St is given 
by: 


Tobs'^t j-vp, . 
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where Toba is the lifetime of the satellite and Vt is 
the relative velocity of source with respect to the lens 
(|Dominik fc Sahu 2000f) . Here, we assume that the as¬ 
trometric centroid shifts more than St = 5-\/2cra are mea¬ 
surable with Gaia in which Ua is the Gaia astrometric 
precision. Since Gaia can measure the astrometric dis¬ 
placement only along the scan, we combine Ua with a 
factor of ^/2 to give the G aia two-dimensional accuracy 
([Belokurov fc Evans 200^ . The amount of the threshold 
impact parameter for a typical astrometric microlensing 
event being observed with Gaia with the astrometric pre¬ 
cision Oa = 300 pas is equal to Ua ^ 25 whereas the 
threshold impact parameter for the photometric obser¬ 
vation of a microlensing event is about one. Hence, the 
astrometric cross-section which is proportional to oc 
is much larger than the photometric one. 


2.1. Astrometric optical depth 

According to the definition of the threshold impact 
parameter Ua, during the observational time Tabs the 
effective area around each lens which yields an astro- 
metric shift in the projected source positions more than 
the threshold amount St is A{x,Mi) = 2uaRETobsVt- 
This area is a function of the lens mass Mi and the lens 
distance Di = xDg from the observer. The astromet¬ 
ric optical depth is defined as the cumulative fraction 
of these areas in the lens plane, for all lens distances 
from the observer to th e source star which is given by 
(iDominik fc Sahu 200011 : 

ra=Ds [ dxpd{x) [ dMi^^Y^A{x,Mi), (4) 
Jo Jo 

where f{Mi) is the distribution function of the lens mass 
and Pdix) is the stellar density distribution in the Galac¬ 
tic disk. We use the follow ing distribution for the Galac¬ 
tic disk dRahal et al. 200911 : 

Pd{x) = — exp(- - -) exp( —), (5) 

where R is the radial distance from the Galactic center, 
z is the altitude with respect to the Galactic plane, H = 
0.325 kpc is the height scale, h = 3.5 kpc is the length 
scale of the disk and E is the column density of the disk at 
the Sun position. By inserting the amount of A{x,Mi)^ 
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the optical depth is: 


Ta = 




f{Mi) 

Vm ■ 


( 6 ) 


The optical depth can be used to estimate the number of 
observed events. If Nbg background source stars brighter 
than 20 G magnitudes are observed during the Gaia life¬ 
time Tabs, the number of source stars with the astromet¬ 
ric shift, in their projected positions, more than St is 
given by: 


Na 


'^ohs^bg 
2 tEUa 


(7) 


By considering e as the efficiency for measuring the lens 
mass of these observed events, the number of astrometric 
microlensing events being observed with Gaia in which 
the lens mass can be measured is Ng = Na e. 


3. BINARY ASTROMETRIC MICROLENSING DURING THE 
GAIA MISSION 

The binary fractions of massive stellar populations 
have not explicitly been specified yet. However, this issue 
is very important for studying the gravitational waves. 
Because, the sources of the gravitational waves are binary 
systems composed of white dwarfs, neutron stars and 
black holes e.g. Riles (2013), Cutler and Thorne (2002). 
Here, we want to investigate whether Gaia can specify 
the binary fractions for these stellar populations through 
binary astrometric microlensing. For this goal, we esti¬ 
mate the number of detectable binary events which have 
considerable binary signals during the Gaia era by per¬ 
forming Monte Carlo simulation according to the Gaia 
observing strategy. 

In the first step, we produce an ensemble of binary 
astrometric and photometric microlensing events due 
to massive stellar populations as microlenses according 
to their physical distributions explained in subsection 
m- We consider binary stellar-mass black holes, 
neutron stars, white dwarfs and main-sequence stars 
as microlenses. Then, corresponding light curves and 
astrometric trajectories are re-generated according to 
the Gaia observing strategy. By c onsi dering a suitable 
criterion explained in subsection (IT^ . we investigate 
whether the binary signatures can be inferred from the 
Gaia observations. Finally, by calculating the optical 
depth (explained in subsection 12.11) . the number of 
detectable binary events during the Gaia mission are 
estimated. The results are gathered in subsection (13.2p . 


3.1. parameter space 

The distribution functions used to generate the binary 
astrometric and photometric microlensing events are ex¬ 
plained here. 

For the source stars, we make an ensemble of the 
source stars according to their mass, age, metallicity and 
magnitude and then choose uniformly the source star 
from that. We take the ages of the stars in the range 
of log{t/yr) = [6.6,10.2] with intervals of Alog(t) = 
0.05 dex. For each age range, we estimate the num¬ 


ber of stars according to a pos i tive star form ation rate 
dN/dt > 0 ([Cignoni et al. 20061 : iKarami 20ldll . For each 
star in this age range, we take the mass of source star, 
i.e. M*, from the Kroupa mass function in the range 
of M* G [O.1,3]M0 (jKroupa et al. 19^ iKroupa 20011) 
and metallicity, i.e. Z, according to the age—metallicity 
relation, and the distribution function for the metallic¬ 
ity from Twarog (1980). Finally, we use the theoreti¬ 
cal isochrones of the colour—magnitude diagram (CMD) 
obtained using the numerical simu lation of the stellar 
evolution from the Padova model (|Marigo et al. 200^ 
to simulate the absolute magnitude of each star. The 
isochrones are defined for stars with different masses, 
colours and magnitudes but with the same metallicities 
and ages. 

The coordinate of the source star toward the Galactic 
bulge (l,b,Ds), is selected from dN/dfl oc p{l,b,Ds)Dl 
in the range of I G [—10,10]°, b G [—5,5]° and 
Ds G [0,10] kpc where p{l,b,Ds) is the combination 
of the stellar dens ities in the Galactic bulge and disk 
(|Rahal et al. 20091) . Galactic bulge has a triaxial struc¬ 
ture rotated along the main axes of bulge according to 
the orientation angles: rj = 78.9° represents the angle 
between the bulge major axis and the line perpendicu¬ 
lar to the Sun-Galactic center line, /3 = 3.5° is the an¬ 
gle between the bulge plane and the Galactic plane and 

J ' = 91.3° is the rol l angle around the major axis of bulge 
Robin et al. 200,^ . 

We obtain the apparent magnitude of the source star 
(m*) according to the absolute magnitude, distance mod¬ 
ulus and extinction. The extinction map towards the 
Galactic bulge and in the Galactic distance is obtained 
in some references, e.g. Gonzalez et al. (2012), Nataf et 
al. (2013). But the extinction map in the closer distances 
is not clear. We know that the extinction can be deter¬ 
mined according to t he column density of Hydrog en con¬ 
tent of the Galaxy ([Weingartner fc Draine 200lh . The 
different fractions of Inter Stellar Medium (ISM) in dif¬ 
ferent directions toward the Galactic bulge are due to Hy¬ 
drogen content. Hence, we need to estimate these frac¬ 
tions in different directions toward the Galactic bulge. 
We do this by comparing (in the each direction at the 
Galactic distance) (i) the extinction map provided by 
Gonzalez et al. (2012) and (ii) the extinction amount 
calculated according to its relation with the Hydrogen 
column density by assuming all ISMs are made from the 
Hydrogen content. We use ISM density law introduced 
by Robin et al. (2003). Also we obtain the V-band ex¬ 
tinction from its value in near infrared bands according 
to the relations introduced by Cardelli et al. (1989). In¬ 
deed, we assume the fraction of ISM composed from the 
Hydrogen content in the Galactic scale depends only on 
the direction but not on the distance. The source stars 
fainter than 20 G magnitudes are not considered. 

For generating the light curve of microlensing events 
with the finite source size effect, we need the radius of 
source star. For the main-sequence stars, the relation 
between the mass and radius is given by 
where all parameters normalized to the sun’s value. We 
do not consider giant source stars. The Gaia astrometric 
and photometric precisions are poin ted out according t o 
the apparent magnitude of source (|Varadi et al. 200^ . 
We consider the events in which 9e > 5-\/2cto. 
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Figure 1. A typical binary microlensing event with binary stellar- 
mass black holes. The light curves (left panel), astrometric trajec¬ 
tories (right panel) and the source trajectories with respect to the 
caustic curve (the inset in the left-hand panel) without and with 
considering the parallax effect are shown with green dashed and 
blue dotted lines respectively. The point-mass lens models (consid¬ 
ering the parallax effect) are shown by red solid lines. Data points 
taken by Gaia are shown with black points. The relevant parame¬ 
ters of the source and lens are M* = 1.1 Mq, m* = 15.8 mag, 
Ds = 2.0 kpc, ^ = —20°, vt = 120 kms~^, uq = 0.4 and 
Ml = 10 Mq, q = 0.6, Di = 0.4 kpc, d = 0.9 Re respectively. 

For indicating the trajectory of the source with respect 
to x-axis, we identify this path with an impact parameter 
and orientation defined by an angle between the trajec¬ 
tory of source star and x-axis i.e. We let this angle 
change in the range of [0,27r]. The impact parameter 
is taken uniformly in the range of uq G {—Ua,Uo\ and 
the corresponding time for uq is chosen uniformly in the 
range of <o S [0, 5] years which is the lifetime of Gaia. 

To simulate the position of the lens in the Galactic 
disk, by considering the same Galactic latitude b and lon¬ 
gitude I as the source, we take its distance with respect to 
the observer from the probability function of astromet¬ 
ric microlen sing detection dT/dx oc pd{x)( \ — xY , where 
X = Di/Ds (iDominik fc Sahu 20001) . 

The mass of the stellar-mass black holes as microlenses 
is taken from the distribution function introduced by 
Farr et al. (2011). We take the mass of the neutron stars 
from a Gaussian function, with an average and disper¬ 
sion of Mq = 1.28 Mq and a = 0.28 Mq in the range of 
[1,2.5] Mq (lOzel et al. 2012I) . The mass of white dwarfs 
is taken from a combination of two Gaussian functions 
with the averages and dispersions of Mq = 0.59 Mq, 
a = 0.047 Mq and Mq = 0.711 Mq, a = 0.109 Mq 
for two different samples of white dwarfs in the range 
of [0.3,1.4] Mq (IKeoler et al. 200^ . We take the mass 
of main-sequence stars as microlenses from the Kroupa 
mass function in the range of [0.1, 2] Mq. 

The velocities of the lens and the source star are taken 
from the combination of the global an d dispersion veloc¬ 
ities of the Galactic disk and bulge (IRahal et al. 20091 : 
Binnev fc Tremaine 1987] ). 

To simulate the parallax effect, we first convert the 
Galactic coordinate of the source to the ecliptic coordi¬ 
nate using LAMBDA sit43- Then, the Earth’s position 
vector with respect to the Sun is calculated (e.g. the Ap¬ 
pendix section from An et al. (2002)). The time when 
the Earth is in the perihelion is = — 18days with re- 

® http://lambda.gsfc.nasa.gov/toolbox/tb-coordconv.cfm 


spect to the time when Gaia entered its operational orbit 
i.e. start time of observation. 

Here, we review the Gaia scanning law. Gaia has two 
astrometric fields of view separated by 106.5°. Its spin 
axis makes an angle of 45° with respect to the direc¬ 
tion towards the Sun. The satellite spins around this 
axis with the angular rate some of 60 as s~^. Hence, 
observations occur in pair whose time interval is 106.5 
minutes. The spin axis in turn has a precession motion 
around a circle. Finally, the spacecraft will move in a 
Lissajous orbit around the second Lagrange point of the 
Sun-Earth system. The combination of these motions 
gives the Gaia scanning law 0. Each object will be ob¬ 
served 72 times averagely with the cadence some of 25 
days. This satellite measures the positions of objects just 
set along the scan direction. Eor simulating the obser¬ 
vational astrometric trajectories and photometric light 
curves, we adopt the Gaussian distribution for the ca¬ 
dence between two pairs of data taken by Gaia with the 
average and dispersion of 25 and 5 days. Each pair oc¬ 
curs in 106.5 minutes. In each observation we consider a 
random scanning direction chosen uniformly in the range 
of [0,27r]. 

Here, we explain distribution functions used for binary 
systems. For all binary systems we take the mass ratio 
from the distribution function introduced by Duquennoy 
& Mayor (1991). We consider binary stellar-mass 
black holes (BBH), binary neutron stars (BNS), binary 
white dwarfs (BWD) and binary main-sequence stars 
(BMS) as microlenses. (Note that, BNS and BWD are 
not necessarily double neutron stars and double white 
dwarfs. Their primary lenses are a neutron star and a 
white dwarf respectively.) We choose the semi-major 
axis s of the binary orbit from the Opik’s law where the 
distribution function for the primary-secondary distance 
is proportional to p{s) = dN/ds oc s~^ (jOpik 19241) 
in the range of s G [0.6,30] A.U.. We use the gen- 
eralized version of the adaptive contouring algorithm 
(IDominik 200711 for simulating astrometric trajectories 
in binary microlensing events. 

We let the microlenses rotate around their center of 
mass. Indeed, the perturbations developed with time 
such as lens orbital motion can make considerable devi¬ 
ations in astrometric t rajectories whic h helps to discover 
the second microlens (jSaiadian 20141) . The parameters 
of lenses orbit are taken as follows: the time of arriving 
at the perihelion point of orbit tp is chosen uniformly in 
the range of tp € [to — + ^*1 where P is the orbital 

period of lenses motion. The projection angles to spec¬ 
ify the orientation of orbit of lenses with respect to the 
sky plane, i.e. S and d, are taken uniformly in the range 
of [—f, f]- We take the eccentricity of the lenses’ orbit 
uniformly in the range of e G [0,1]. 

3.2. results 

To determine the detectability of binary signals we 
calculate x^- Two models are fitted to the astromet¬ 
ric and photometric data: (a) the known rotating bi¬ 
nary microlensing model with the parameters used to 
simulate data, i.e. xl (b) the simple Paczyhski mi- 

^ http://sci.esa.int/gaia/40846-galactic-structure 
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Table 1 

The results of the simulation of the binary astrometric microlensing events with GAIA. 



tsidays) 

tae(days) 

Ua 

eb{%) 

S {Mq pc-^) 

7a X 10-’’ 

Nb{f = 0.1) 

Nbif = 0.3) 

Nbif = 0.5) 

BBH 

90.2 

568.6 

18.3 

13.5 

2.09 

8.0 

5.6 

16.9 

28.1 

BNS 

27.1 

109.5 

32.8 

16.5 

0.91 

8.6 

11.4 

34.2 

57.1 

BWD 

17.9 

66.0 

34.1 

15.6 

2.41 

34.8 

76.8 

230.3 

383.8 

BMS 

11.5 

42.5 

43.3 

17.6 

21.3 

434.9 

1316.1 

3948.4 

6580.6 


Note. — The average values of the Einstein crossing time tE{days), the duration of the astrometric event tae (days) threshold impact 
parameter UT, efficiency of detecting binary signatures of microlenses €b{%), the column density E {Mq pc~^), astrometrical optical depth 
tT and number of binary astrometric microlensing events detected with Gaia due to different binary systems: binary stellar-mass black 
holes (BBH), binary neutron stars (BNS), binary white dwarfs (BWD) and binary main-sequence stars (BMS) considering different binary 
fractions /. 


crolensing events with the same Paczyhski parameters, 
i.e. Xp- Our criterion for detectability of the binary sig- 
nal is Ax^ = Xp - xl > ^Xth- Xj (for j=b,P) from 
fitting the astrometric trajectory along the scan Og and 
the magnification factor A is given by: 


2 _ 

Xj — ^i=l I 2 


{0s,i - 


( 8 ) 


where N is the number of data points and cta is the 
Gaia photometric precision. We assume that from fit¬ 
ting process and searching all parameter space the best- 
fitted solution is the known binary solution. We consider 

Ax?. = 200. 

In Figure m we show a typical simulated binary mi¬ 
crolensing event with binary stellar-mass black holes be¬ 
ing observed with Gaia. The light curves (left panel), 
astrometric trajectories (right panel) and the source tra¬ 
jectories with respect to the caustic curve (the inset in 
the left-hand panel) without and with considering the 
parallax effect are shown with green dashed and blue 
dotted lines respectively. The point-mass lens models 
are shown by red solid lines. Data points taken by Gaia 
are shown with black points. The relevant parameters of 
the source and lens are M* = 1.1 Mq, nrii, = 15.8 mag, 
Dg = 2.0 kpc, ^ = —20°, Vt = 120 kms~^, uq = 0.4 and 
Ml = 10 Mq, q = 0.6, Di = 0.4 kpc, d — 0.9 Re respec¬ 
tively. In this event, the binary signal is detectable with 
Gaia. Although, in the photometric light curve there 
is one data point with significant deviation with respect 
to the simple microlensing light curves but in the astro- 
metric trajectory there are several deviated data points. 
Hence, although the Gaia cadence is too sparse so that 
some of caustic-crossing features may be missed in the 
photometric light curves but the astrometric measure¬ 
ments can much likely resolve the binary signatures from 
the second microlenses. Indeed, the astrometric cross 
section is much larger than the photometric one. 

In Table (P) we provide the average values of the Ein¬ 
stein crossing time ts, the duration of the astrometric 
event tae, threshold impact parameter u^, efficiency for 
detecting binary signals e., i.e. the fraction of the sim¬ 
ulated binary events in which the binary signals are de¬ 
tectable with Ax^ > ^Xth given in per cent, the col¬ 
umn densities for different stellar populations E (taken 
from Flynn et al. (2006) and Chabrier (2001)) and the 
astrometrical optical depth from the second to sixth 


column due to different binary systems. Note that, we 
calculate the optical depth for different directions toward 
the Galactic bulge, i.e. I in the range of [—10,10]° and b 
in the range of [—5,5]° and average over all of them. 

The expected number of binary microlensing events 
being observed with Gaia with detectable binary signals, 
i.e. Nb, is given by: 

Nb = Na f Eb, (9) 

where / is the fraction of binary objects in different stel¬ 
lar populations. To calculate Na (equation [7]) we insert 
the average values of Ie, ^ for each stellar population. 
Also, we set the number of background s ource stars to - 
wards the Galactic bulge Nbg = 3.0 x 10® (|Caciari 201411 . 
Noting that the number of background stars brighter 
than 20-G magnitudes for whole sky is on the order of 
billion objects. The amounts of Nb for three values of 
/ = 0.1, 0.3 and 0.5 are reported in Table O- 

According to this table: (i) The Gaia efficiency for de¬ 
tecting different binary systems as microlenses is almost 
10 — 20 per cent which is acceptable according to the 
sparse cadence of Gaia. It means that Gaia detects the 
signal from the second lens with probability ^ 10 — 20 
per cent. 

(ii) The Gaia efficiency decreases with increasing the 
lens mass. Because, by increasing the lens mass, the 
normalized distance between two lens decreases. Hence, 
binary black holes usually make close binaries for which 
the probability of caustic crossing is smaller than that 
in the intermediate binaries due to the less massive mi¬ 
crolenses. 

(iii) Gaia can potentially determine the binary fraction 
of massive stellar populations. Because, the numbers of 
estimated events for different binary fractions are larger 
than one. The number of estimated events for binary 
black holes is small which means that the uncertainty 
for specifying their binary fraction is high, specially if 
the real binary fraction would be less than 0.1. 

(i v) The duration of the astrometric events is give n 
by ([Dominik fc Sahu 20001 : iBelokurov fc Evans 200211 : 
tae = tE 6 e/{^V^ cTa) which is proportional to the Ein¬ 
stein crossing time. Indeed, the ratio of 6El5V2aa is al¬ 
most constant for different stellar populations and larger 
than one. Hence, the more massive microlenses with 
longer Einstein crossing time have the longer astromet¬ 
ric durations. The longest duration of the astrometric 
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events is that due to binary stellar-mass black holes and 
some of 1.56 yrs. This time is one third of the Gaia 
lifetime. The more massive microlenses have longer as¬ 
trometric duration, comparable with the Gaia lifetime 
which are not suitable to be observed with Gaia. 

4. MEASURING THE MASS OF MASSIVE STELLAR 
POPULATIONS 

Here we investigate how Gaia is efficient to measure 
the mass of the massive stellar populations through as¬ 
trometric microlensing of single lenses and specify their 
mass distributions. In this regard, we quantitatively esti¬ 
mate the Gaia efficiency and the number of high-quality 
astrometric microlensing events by performing a Monte 
Carlo simulation. 

We consider different stellar populations, i.e. stellar- 
mass black holes, neutron stars, white dwarfs and 
main-sequence stars as microlenses and simulate as¬ 
trometric microlensing events due to these microlenses 
according to their distribution functions explained in 
subsection (EH). Having re-generated the astrometric 
and photometric curves with synthetic data points 
according to the Gaia observing strategy, we should 
indicate whether the lens mass can be inferred from 
data. In that case, we evaluate the uncertainty in 
the lens mass i.e. ctm, in the same way as Belokurov 
& Evans (2002), and accept the events with the rela¬ 
tive error i.e. um/M i less than 0.5 as high-quality events. 

4.1. results 

The results of this Monte Carlo simulation are summa¬ 
rized in Table ([2|) which lists the average values of the 
lens mass Mi (Mq), Einstein crossing time tE {days), 
duration of the astrometric event tae (days), thresh¬ 
old impact parameter u^, lens-source transverse velocity 
Vt {km s~^), source distance from the observer {kpc), 
astrometric precision of Gaia {y, as), efficiency for de¬ 
tecting high-quality events e (%), i.e. the number of the 
high-quality events divided by the total number of simu¬ 
lated events given in per cent and finally the number of 
high-quality events W (explained in subsection 12.11) for 
four different stellar populations. 

Interpreting this table: (i) The more massive lenses, 
the larger efficiencies of measuring the lens mass. In¬ 
deed, the Einstein crossing time in astrometric microlens¬ 
ing events due to massive lenses is long. As a result, the 
number of data over the astrometric trajectories and pho¬ 
tometric light curves increases. On the other hand, their 
angular Einstein radius and astrometric signal are high. 

(ii) The Gaia astrometric precision for astrometric mi¬ 
crolensing due to more massive microlenses is larger than 
that due to less massive ones. Indeed, we consider two re¬ 
strictions for simulating astrometric microlensing events: 
(a) the source star should be brighter than 20 G magni¬ 
tude and (b) the angular Einstein radius should be larger 
than 5-\/2CTa. Hence, the average value of CTo, reported 
in seventh column, represents the minimum amount of 
the angular Einstein radius which in turn is an increas¬ 
ing function with respect to the lens mass. As a result, 
the threshold impact parameter which is proportional to 
I/'/St decreases with enhancing the lens mass. 

(iii) The number of high-quality events due to the 
massive stellar populations are larger than one. Hence, 


we can obtain some information about their masses and 
mass distribution. A significant point is that this method 
can even give the mass of isolated objects. For the main- 
sequence stars as microlenses, we conclude that some of 
786 high-quality events toward the Galactic bulge will 
probably be confirmed during the Gaia era. However, 
Belokurov and Evans (2002) estimated that number of 
high-quality events due to main-sequence stars for whole 
sky is ^ 2500. Indeed, we consider the number of back¬ 
ground stars brighter than 20 G magnitudes towards the 
Galactic bulge as 300 million objects which is so fewer 
than that was considered by them, i.e. one billion ob¬ 
jects. Meantime, there are some differences between the 
distribution functions used in our simulation and those 
used in the referred paper. However, they estimated just 
high-quality events due to main-sequence stars. 

Note that the mass ranges of these four stellar popula¬ 
tions as microlenses overlap in edges. But, we can discern 
the type of microlenses with the masses in the common 
ranges (belong to two populations) through their color 
and magnitude. Because, these microlenses most likely 
locate at the distance less than I kilo parsec from the 
observer. 

In order to study the dependance of efficiency on the 
parameters of the model, in Figure (H we plot the nor¬ 
malized efficiencies, i.e. e„, for detecting high-quality 
events due to stellar-mass black holes (black solid line), 
neutron stars (red dashed line), white dwarfs (green dot¬ 
ted line) and main-sequence stars (blue dot-dashed line) 
in terms of the relevant parameters of the lens and source. 
To compare the efficiency curves with each other, we nor¬ 
malize all efficiency curves to one which means the area 
under each curve is one, hence we name them as nor¬ 
malized efficiency curves. Since the area under each effi¬ 
ciency curve is proportional to the related value reported 
in the ninth column of Table ([2]) i.e. e, the normalization 
is done by dividing Y values to e. Generally, the behav¬ 
iors of efficiency for different stellar populations are simi¬ 
lar, although some differences exist in the slope of curves, 
peak position, etc. The normalized efficiency function in 
terms of the lens and source parameters is given as fol¬ 
lows: 

(i) The first one is the impact parameter, uq. The 
efficiency of detecting high-quality events is high for as¬ 
trometric microlensing events with uq < 1 for which the 
photometric microlensing can be observed as well as the 
astrometric one. 

(ii) The second parameter is the relative velocity of 
source with respect to lens, Vt- The efficiency increases 
with decreasing the relative velocity. Indeed, the Ein¬ 
stein crossing time rises with decreasing the relative ve¬ 
locity. Since the Gaia cadence is almost constant for 
all events towards the Galactic bulge, so the number of 
observational data for long-duration microlensing events 
(over the astrometric trajectories and photometric light 
curves) is high. As a result, the errors in parameters for 
these events are smaller than those due to short-duration 
ones. 

(iii) and (iv) Two next parameters: M* the mass and 
m* apparent magnitude of source star. The Gaia astro- 
metric precision depends strongly on the source magni¬ 
tude and increases with decreasing it. Hence, the Gaia 
efficiency for measuring the lens mass in astrometric mi¬ 
crolensing due to brighter and heavier sources is higher 
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Table 2 

The results of the simulation of astrometric microlensing events due to massive stellar populations with GAIA. 



Ml (Mq) 

ts (days) 

tae (days) 

Ua 

Vt {km s 

Ds (kpc) 

(Ta (y as) 

£(%) 

Ne 

(BH) 

10.3 

74.1 

460.2 

20.1 

179.8 

4.3 

478.8 

9.8 

44.5 

(NS) 

1.4 

20.5 

86.7 

29.0 

182.5 

3.5 

324.2 

2.9 

34.2 

(WD) 

0.6 

12.3 

46.6 

35.0 

182.7 

3.1 

268.9 

1.2 

76.3 

(MS) 

0.4 

9.1 

33.3 

41.0 

182.8 

2.7 

232.1 

0.8 

785.8 


Note. — The average values of the lens mass M; (Mq), Einstein crossing time ts (days), duration of the astrometric event tae (days), 
threshold impact parameter u^, lens-source transverse velocity vt (km source distance from the observer Ds (kpc), astrometric 

precision of Gaia oT (fias), detection efficiency for high-quality events e (%) and finally the number of high-quality events being observed 
during the Gaia era for four different stellar populations: stellar-mass black holes (BH), neutron stars (NS), white dwarfs (WD) and 
main-sequence stars (MS). 



Figure 2. The normalized efficiency curves, of Gaia for measuring the lens mass in astrometric microlensing due to stellar-mass black 
holes (black solid line), neutron stars (red dashed line), white dwarfs (green dotted line) and main-sequence stars (blue dot-dashed line) 
for some parameters of source star and lens. Noting that all efficiency curves are normalized to one which means the area under each curve 
is one, so-called as normalized efficiency curves. Also, the lens mass for each stellar population is divided to its maximum amount so that 
the range of the normalized lens mass for all populations is in the range of [0,1]. 


than that due to fainter and less massive ones. 

(v) The next parameter is the normalized lens mass, 
Ml. Since the mass range for all stellar population is not 
similar, we normalize the lens mass for each stellar pop¬ 
ulation to its maximum amount so that the range of the 
normalized lens mass for all populations is in the range of 
[0,1]. For the massive lenses the Einstein crossing time 
is long. Consequently, the number of data over the astro- 
metric trajectory and photometric light curve increases. 
On the other hand, their angular Einstein radius and as¬ 
trometric signal are high. Therefore, the Gaia efficiency 
for measuring the mass of more massive lenses is higher 
than that of less massive lenses. 

(vi) The next parameter is the corresponding time to 
the closest approach, to- We uniformly choose to in the 
range of [0, 5]T. But the events with to 2.5T has the 
maximum efficiency. Because, for these events the num¬ 


ber of observational data over the astrometric trajectory 
and light curve is maximum. 

(vii) and (viii) Two last parameters: x = Di/Dg the 
relative distance of the lens to the source star from the 
observer and Dg the source distance from observer. In¬ 
deed there are two kinds of suitable sources for the Gaia 
detection: (a) nearby sources with Dg < \ kpc and (b) 
the bright sources often located at the Galactic disk. 
Correspondingly, there are two peaks in the efficiency 
diagram versus Dg. These two classes have the different 
xs. We know that the Einstein crossing time is propor¬ 
tional to is oc Dgx{l — x), but the relative parallax 
and the angular Einstein radius depend on these param¬ 
eters as 9E oa oc ^ Although, the efficiency 

enhances with increasing each of them, but the first is an 
increasing function versus x and Dg and the others are 
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Figure 3. The normalized efficiency of Gaia for measuring the 
lens mass in astrometric microlensing due to different stellar pop¬ 
ulations (top panel) and the average value of the extinction in V- 
band, Av, (red solid line) versus the Galactic latitude b. Note that 
the efficiency curve for each population normalized to the related 
cumulative value reported in the ninth column of Table [2] 

decreasing functions. Therefore, the efficiency is maxi- 
mized for the first class when x ^ [0.6,0.7] and for the 
second class when x <C 1. 

According to the dependance of efficiency on the pa¬ 
rameters, we expect that most high-quality events have 
long Einstein crossing times or small impact parameters 
or their source is a bright star. Hence, most sources are 
so close to the observer. 

We also investigate the dependence of the Gaia (nor¬ 
malized) efficiency on the Galactic latitude, shown in the 
top panel of Figure ([31). The normalized efficiency de¬ 
creases very slowly by increasing the distance from the 
Galactic plane. Because, there are two effects act re¬ 
versely: (i) The extinction toward the Galactic bulge 
strongly increases by decreasing the distance from the 
Galactic plane as shown in the bottom panel of Figure ([3]) 
and (ii) the Galactic disk density, equation decreases 
by increasing distance from the Galactic plane. The first 
effect decreases the probability of finding a source star 
brighter than 20-G magnitude toward the Galactic cen¬ 
ter while the abundance of stars is maximum in this di¬ 
rection. Gonsequently, the efficiency is almost constant 
over the Galactic bulge. To determine the extinction, we 
average over the amounts of the extinction due to the 
Galactic longitudes in the range of [—10,10]° for each 
value of b. 

5. CONCLUSIONS 

One of astronomical phenomena being observed with 
Gaia is astrometric microlensing events. In this work, we 
investigated whether Gaia can specify the binary frac¬ 
tions and mass distributions of massive stellar popula¬ 
tions through astrometric microlensing. 

We simulated binary astrometric microlensing events 
being observed with Gaia due to binary stellar-mass 
black holes, binary neutron stars, binary white dwarfs 
and binary main-sequence stars to evaluate the effi¬ 
ciency for detecting the binary signatures. We consid¬ 
ered > 200 as the detectability criterion and con¬ 
cluded that Gaia detects the signal from the second lens 
with probability ~ 10 — 20 per cent. Then, We esti¬ 
mated the number of binary astrometric microlensing 


due to the mentioned populations with detectable binary 
signals in the Gaia era as 6, 11, 77 and 1316 respec¬ 
tively (/ = 0.1). Hence, Gaia can potentially determine 
the binary fraction of massive stellar populations. Since 
the predicted number of binary astrometric microlensing 
events of stellar-mass black holes is small, so the Gaia 
accuracy to specify their binary fraction is low. The sig¬ 
nificant points is that binary systems composed of white 
dwarfs, neutron stars or black holes are sources of de¬ 
tectable gravitational waves. Therefore, indicating the 
binary fraction of these massive stellar populations helps 
to assess the amount of the gravitational waves. 

Also, we investigated the Gaia efficiency for measur¬ 
ing the mass of massive stellar populations in the Galac¬ 
tic disk, i.e. stellar-mass black holes, neutron stars and 
white dwarfs in addition to main-sequence stars as mi¬ 
crolenses in astrometric microlensing events of single 
lenses. The Gaia efficiency for measuring the mass of 
microlenses enhances with increasing the lens mass. By 
performing Monte Garlo simulation, we concluded that 
the Gaia efficiencies for detecting high-quality events due 
to these stellar populations are 9.8, 2.9, 1.2 and 0.8 per 
cent respectively. We estimated the number of high- 
quality events for these populations as 45, 34, 76 and 
786 respectively. Hence, Gaia can obtain some informa¬ 
tion about their masses and mass distribution of these 
massive stellar populations. 

I am grateful to Sohrab Rahvar for helpful discussions 
and comments and Scott Gaudi for careful reading and 
commenting on the manuscript. I thank anonymous ref¬ 
eree for useful comments and suggestions which certainly 
improved the manuscript. Also, I acknowledge the use 
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